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ABSTRACT 

We have observed the Herbig Ae/Be sources AB Aur, VV Ser, V1685 Cyg 
(BD+40°4124), AS 442, and MWC 1080 with the Palomar Testbed Interferom- 
eter, obtaining the longest baseline near-IR interferometric observations of this 
class of objects. All of the sources are resolved at 2.2 /zm with angular size 
scales generally < 5 mas, consistent with the only previous near-IR interfero- 
metric measurements of Herbig Ae/Be stars by Millan-Gabet and collaborators. 
We determine the angular size scales and orientations predicted by uniform disk, 
Gaussian, ring, and accretion disk models. Although it is difficult to distinguish 
different radial distributions, we are able to place firm constraints on the incli- 
nations of these models, and our measurements are the first that show evidence 
for significantly inclined morphologies. In addition, the derived angular sizes for 
the early type Herbig Be stars in our sample, V1685 Cyg and MWC 1080, agree 
reasonably well with those predicted by the face-on accretion disk models used by 
Hillenbrand and collaborators to explain observed spectral energy distributions. 
In contrast, our data for the later- type sources AB Aur, VV Ser, and AS 442 are 
somewhat inconsistent with these models, and may be explained better through 
the puffed-up inner disk models of Dullemond and collaborators. 
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1. Introduction 

Herbig Ae/Be (HAEBE; Herbig 1960) stars are intermediate- mass (2-10 M Q ) young 
stellar objects that show broad emission lines, rapid variability, and excess infrared and 
millimeter-wavelength emission. These properties are consistent with the presence of hot 
and cold circumstellar dust and gas. While there is still some debate about the morphology 
of this circumstellar material, most evidence supports the hypothesis that in many cases the 
dust and gas lies in a massive (~ 0.01 M ) circumstellar disk (Natta, Grinin, & Mannings 
2000; Hillenbrand et al. 1992, hereafter HSVK). 

The strongest evidence for circumstellar disks around HAEBE stars comes from direct 
imaging with millimeter interferometry. Flattened structures around several sources have 
been resolved on ~ 100 AU scales (Mannings & Sargent 1997; Mannings & Sargent 2000; 
Pietu et al. 2003), and detailed kinematic modeling of one source, MWC 480, shows that 
the observations are fit well by a rotating Keplerian disk (Mannings et al. 1997). For a 
spherical distribution, these and other observations (e.g., Mannings 1994) imply extinctions 
at visible and infra-red wavelengths much higher than actually observed. In addition, in 
recent Ha spectropolarimetric observations of HAEBE sources (which trace dust on scales of 
tens of stellar radii), Vink et al. (2003) find signatures of flattened circumstellar structures 
around 83% of their sample, and evidence for rotation around 9 HAe stars. Furthermore, the 
forbidden emission lines that arise in winds and outflows around HAEBE sources typically 
show blue-shifted emission but lack redshifted emission, which suggests that the redshifted 
component of the outflow is occluded by a circumstellar disk. The broad linewidths of low- 
velocity features are consistent with this emission arising in rotating circumstellar disk winds 
(Corcoran & Ray 1997). 

The distribution of circumstellar material around HAEBEs can also be inferred from 
modeling of spectral energy distributions (SEDs). Three distinct morphologies were iden- 
tified in this way by HSVK, who classified observed HAEBE sources into three groups, I, 
II, and III. All sources in our observed sample fall into Group I, which has SEDs of the 
form \F\ oc A 4 / 3 . These can be modeled well by flat, irradiated, accretion disks with in- 
ner holes on the order of ~ 10 stellar radii. Recent SED modeling of a sample of fourteen 
isolated HAEBE stars with the characteristics of Group I sources is consistent with emis- 
sion from a passive reprocessing disk (Meeus et al. 2001). Moreover, Meeus et al. (2001) 
(and other investigators, e.g., Natta et al. 2001) attribute this emission to the outer part 
of a flared circumstellar disk (e.g., Chiang & Goldreich 1997), while previous authors at- 
tributed blackbody components observed in SEDs of HAEBE sources to tenuous envelopes 
(Hartmann et al. 1993; Miroshnichenko et al. 1999; Natta et al. 1993). 

Size scales and orientations of disks around HAEBE stars can only be determined di- 
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rectly through high angular resolution imaging. The spatial and velocity structure of cooler 
outer HAEBE disks on ~ 100 AU scales has been mapped with millimeter-wave interferom- 
eters (as discussed above). To probe the warmer inner regions of the disk (~ 1 AU scales), 
measurements with near-IR interferometers are necessary. The only near-IR interferometric 
observations of HAEBE sources to date, conducted with the IOTA interferometer (Millan- 
Gabet et al. 1999; Millan-Gabet, Schloerb, & Traub 2001, hereafter MST), led to sizes and 
orientations of sources largely inconsistent with values estimated using other techniques. 
A geometrically flat disk may be too simplistic to accommodate all the observations, and 
puffed up inner disk walls (Dullemond, Dominik, & Natta 2001; hereafter DDN) or flared 
outer disks (e.g., Chiang & Goldreich 1997) may need to be included in the models. How- 
ever, the limited number of HAEBE sources observed with near-IR interferometers and the 
sparse u — v coverage of these observations (MST) make it difficult to draw unambiguous 
conclusions about the structure of the circumstellar material. 

We have begun a program with the Palomar Testbed Interferometer (PTI) to observe 
HAEBE stars. By increasing the sample size and improving u — v coverage, we aim to 
understand better the structure of the circumstellar emission on ~ 1 AU scales. In this 
paper, we present results for five sources, AB Aur, VV Ser, V1685 Cyg (BD+40°4124), AS 
442, and MWC 1080. We note that neither our list of HAEBEs, nor that of MST, represents 
an unbiased sample, but rather, is limited to those stars that are bright enough (K < 6.5) to 
be successfully observed. We model the structure of circumstellar dust around HAEBE stars 
using these PTI data, together with IOTA measurements where available. Specifically, we 
compare various models — Gaussians, uniform disks, uniform rings, and accretion disks with 
inner holes — to the visibility data to determine approximate size scales and orientations of 
the circumstellar emission. 

In §2, we describe the PTI observations. In §3, we fit the observed data to several differ- 
ent models for the circumstellar dust distribution and derive angular sizes and orientations. 
Implications of the modeling and comparisons with previous observations are discussed in 
§4. 

2. Observations and Calibration 

The Palomar Testbed Interferometer (PTI) is a long-baseline near-IR Michelson inter- 
ferometer located on Palomar Mountain near San Diego, CA (Colavita et al. 1999). PTI 
combines starlight from two 40-cm aperture telescopes using a Michelson beam combiner, 
and records the resulting fringe visibilities. These fringe visibilities are related to the source 
brightness distribution via the van Cittert-Zernike theorem, which states that the visibility 
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distribution in u — v space and the brightness distribution on the sky are Fourier transform 
pairs (Born & Wolf 1999). 

We observed five HAEBE sources, AB Aur, VV Ser, V1685 Cyg (BD+40°4124), AS 442, 
and MWC 1080, with PTI between May and October of 2002. Properties of the sample are 
included in Table 1. We obtained K-band (2.2 /iin) measurements on an 85-m North-West 
(NW) baseline for all five objects, and on a 110-m North-South (NS) baseline for three. The 
NW baseline is oriented 109° west of north and has a fringe spacing of ~ 5 mas, and the 
NS baseline is 160° west of north and has a fringe spacing of ~ 4 mas. A summary of the 
observations is given in Table 2. 

PTI measures fringes in two channels, corresponding to the two outputs from the beam 
combiner. One output is spatially filtered with an optical fiber and dispersed onto five 
"spectral" pixels, while the other output is focused onto a single "white-light" pixel (without 
spatial filtering). The white-light pixel is used principally for fringe-tracking: the fringe phase 
is measured and then used to control the delay line system to track atmospheric fringe motion 
(and thus maintain zero optical path difference between the two interfering beams). The 
spectral pixels are generally used to make accurate measurements of the squared visibility 
amplitudes (V 2 ) of observed sources. We sample the data at either 20 or 50 milliseconds in 
order to make measurements on a timescale shorter than the atmospheric coherence time. 
A "scan", which is the unit of data we will use in the analysis below, consists of 130 seconds 
of data, divided into five equal time blocks. The V 2 is calculated for each of these blocks 
using an incoherent average of the constituent 20 or 50-ms measurements from a synthetic 
wide-band channel formed from the five spectral pixels (Colavita 1999). V 2 for the entire 
scan is given by the mean of these five estimates, and the statistical uncertainty is given by 
the standard deviation from the mean value. 

We calibrate the measured V 2 for the observed HAEBE sources by comparing them to 
visibilities measured for calibrator sources of known angular sizes, for which we can easily 
calculate the expected V 2 for an ideal system. The visibilities are normalized such that 
V 2 = 1 for a point source observed with an ideal system. We calculate the expected V 2 by 
assuming that the calibrators are uniform stellar disks. Making use of the van Cittert-Zernike 
theorem, the squared visibilities for these sources are given by 

Here, J\ is the first-order Bessel function. 9 is the angular diameter of the star, and r uv is 
the "uv radius" , defined by 

r uv = V^T^=^, (2) 



- 5 - 



where B is the baseline vector, s is a unit vector pointing from the center of the baseline 
towards the source, and A is the observing wavelength. (The qualitative explanation of 
Equation 1 is that while for unresolved sources the visibility is constant with increasing 
uv radius, for progressively larger sources the visibility decreases faster with increasing uv 
radius.) By comparing V^ c to the measured V 2 for a calibrator, we derive the "system 
visibility" , which represents the point source response of the interferometer: 

V 2 

t|2 _ meas, calibrator /q\ 
^sys = 772 • \ 6 > 

calc 

This system visibility, in turn, is used to calibrate the squared visibilities for the target 
source: 

V 2 

tt2 mcas, target / . \ 

^target — 772 ' ^ > 

sys 

Specifically, we determine V^ s at the time of each target scan, using an average of V^ eas ca iib ra tor 
weighted by the proximity of the target and calibrator in both time and angle. For further 
discussion of the calibration procedure, see Boden et al. (1998). 

Calibrators must be close to the target sources (on the sky an in time) so that the 
atmospheric effects will be the same for both. They should also be of small angular size, 0, 
so that V 2 — > 1 and dV 2 ys /d0 — > 0, and the calibration is thus less sensitive to uncertainties 
in the assumed calibrator diameter. The angular size of a calibrator can be estimated from 
the published stellar luminosity and distance, from a blackbody fit to published photometric 
data with the temperature constrained to that expected for the published spectral type, 
or from an unconstrained blackbody fit to the photometric data. We adopt the average of 
these three size estimates in our analysis, and the uncertainty is given by the spread of these 
values. Relevant properties of the calibrators used in these observations are given in Table 
3. 



3. Results 

We measured calibrated squared visibilities for AB Aur, VV Ser, V1685 Cyg, AS 442, 
and MWC 1080 (Table 2). All five sources are resolved by PTI (i.e., V 2 is significantly 
different from unity), implying angular sizes > 1 mas. The data are consistent with disk- 
like morphologies for all sources, and we can place good constraints on disk inclinations for 
most sources. MWC 1080, V1685 Cyg, and VV Ser show evidence for significantly non-zero 
inclinations, while a circularly symmetric distribution appears appropriate for AB Aur. The 
AS 442 data are insufficient to constrain the inclination. 
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Interferometric observations of AB Aur and MWC 1080 at 2.2 /im have also been ob- 
tained with the 21-m and 38-m baselines of the IOTA interferometer (Millan-Gabet et al. 1999; 
MST). When combined with our longer baseline PTI data (85-m and 110-m), these help fill 
in the u — v plane, and enable us to improve constraints on source models (see below; Figure 
11). Based on discussion with R. Millan-Gabet, we assign an uncertainty to each IOTA 
visibility given by the standard deviation of all data obtained for a given source with a given 
baseline. We verify the registration of the IOTA and PTI data using calibrators observed 
by both interferometers. Since the IOTA data for AB Aur was calibrated using HD 32406, 
which is unresolved by both PTI and IOTA, we can be confident of the registration. We 
measured the diameter of HD 220074, the calibrator for MWC 1080, to be # UD = 1.98 ±0.06, 
while MST assumed a size of 2.10 ± 0.22. This difference in angular size translates into only 
a 0.7% effect, which is within the measurement errors (the effect is so small because the 
calibrator is essentially unresolved by IOTA). 



Nearby companions that lie outside the interferometric field of view, ~ 50 mas, but 
within the field of view of the detector, ~ 1", will contribute incoherent light to the visibil- 
ities. For MWC 1080, which has a known nearby companion (Corporon 1998), we use the 
correction factor 



where AK is the difference in K-band magnitudes between the two stars. For MWC 1080, 
we measured AK = 2.70 (angular separation = 0"78) using the Palomar Adaptive Optics 
system on the 200-inch telescope on November 18, 2002. V1685 Cyg is also known to have 
a faint companion (AK = 5.50; Corporon 1998), but the effect of this companion on the 
visibilities is negligible. Adaptive optics images of the other sources in our sample show that 
none of these has any bright companions (AK < 5) at distances between ~ 50 mas and 1". 

Our measured visibilities contain information about emission from both the circumstellar 
material and the star itself. We can remove the effect of the central star on the visibilities 
by including it in the models: 



where is the stellar flux, F x is the excess flux (both measured at 2.2 pm), V* « 1 is the 
visibility of the (unresolved) central star, and V x is the visibility due to the circumstellar 
component. It is reasonable to assume that K ~ 1, since for typical stellar radii (~ 5 R ) 



3.1. Visibility Corrections 




(5) 
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and distances (~ 500 pc), the angular diameters of the central stars will be ~ 0.5 mas. In the 
case of the binary model described below (§3.3.5), we do not perform any such correction, 
since the basic model already includes the stellar component. 

Equation 6 assumes that the central star is a point source, and thus contributes coher- 
ently to the visibilities. It is also possible that the starlight is actually observable only as 
scattered light emission, and that it will have some incoherent contribution to the visibility 
(K 7^ 1). For example, coronographic imaging with the Hubble Space Telescope has re- 
vealed scattered light on angular scales from ~ 0'.'l-9" around AB Aur (Grady et al. 1999). 
A proper treatment of the effects of this scattered light on the visibilities is beyond the scope 
of this work, but we mention it as a possible source of uncertainty. Since the near-IR excess 
from HAEBE sources typically dominates over the near-IR stellar emission (§3.2; Table 1), 
the effect should be insignificant. 



3.2. Photometry 

F* and F x affect the visibilities (Equation 6), and thus it is important to determine these 
quantities accurately. Since HAEBE objects are often highly variable at near-IR wavelengths 
(e.g., Skrutskie et al. 1996), we obtained photometric K-band measurements of the sources 
in our sample that are nearly contemporaneous with our PTI observations, using the Palomar 
200-inch telescope between November 14 and 18, 2002. Calibration relied on observations of 
2MASS sources close in angle to the target sources, and we estimate uncertainties of ~ 0.1 
magnitudes. Our photometry is consistent with published measurements to within ~ 0.3 
magnitudes for all objects (HSVK; Eiroa et al. 2001). 

Following HSVK and MST, we calculate and F x using our K-band photometry (Ta- 
ble 1) combined with BVRI photometry, visual extinctions, and stellar effective temperatures 
from the literature (HSVK; Oudmaijer et al. 2001; Eiroa et al. 2001; Bigay & Gamier 1970). 
De-reddening uses the extinction law of Steenman & The (1991). Assuming that all of the 
short-wavelength flux is due to the central star, we fit a blackbody at the assumed effective 
temperature to the de-reddened BVRI data. The K-band stellar flux is derived from the 
value of this blackbody curve at 2.2 fim. The excess flux is then given by the difference 
between the de-reddened observed flux and the stellar flux. The derived fluxes are given in 
Table 1. We note that VV Ser and AS 442 are optically variable by AV > 1 magnitudes on 
timescales of days to months (while the other sources in our sample show little or no optical 
variability; Herbst & Shevchenko 1999), and thus F, is somewhat uncertain. However, since 
Fx/F* 1 for these objects, this uncertainty is negligible when modeling the visibilities. 
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3.3. Models 

For each source, we compare the observed visibilities to those derived from a uniform 
disk model, a Gaussian model, a ring model, and an accretion disk model with an inner disk 
hole (all models are 2-D). If we assume that the inclination of the circumstellar material is 
zero, then the one remaining free parameter in the models is the angular size scale, 9. When 
we include inclination effects, we fit for three parameters: size (#), inclination angle (0), 
and position angle (ip). Inclination is defined such that a face-on disk has = 0, and ip is 
measured east of north. Following MST, we include <p and ip in our models of the brightness 
distribution via a simple coordinate transformation: 

/ . , . , , ycos^-xsin^ 
x — xsmip + ycosip; y = . (7) 

COS0 

Here, (x, y) are the coordinates on the sky, and (x', y') are the transformed coordinates. 
The effect of this coordinate transformation on the visibilities will be to transform (u, v) to 
(«>'): 

u' — u sin ip + v cos ip; v' — cos (p(v sin ip — u cos ip) . (8) 

Substitution of (x',y') for (x,y), and (u',v') for (u,v) in the expressions below yield models 
with inclination effects included. 

In addition to these four models, we also examine whether the data are consistent with 
a wide binary model, which we approximate with two stationary point sources. For this 
model, the free parameters are the angular separation (9), the position angle (ip), and the 
brightness ratio of the two components (R). 



3.3.1. Gaussian Model 

The brightness distribution for a normalized Gaussian model is given by 

T , \ ( 41n2 (x 2 + y 2 )\ 

Igauss{x, y) = exp I — j , (9) 

and the (normalized) visibilities expected for this observed brightness distribution are cal- 
culated via a Fourier transform to be, 

Wr uv )=exp(-^^). (10) 

Here, (x,y) are the angular offsets from the central star, 9 is the angular FWHM of the 
brightness distribution, and r uv = (u 2 + v 2 ) 1 / 2 is the "uv radius" (Equation 2). The model 
for the observed squared visibilities is obtained by using Equation 6 with V x = Vg auss . 
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3.3.2. Uniform Disk Model 

The brightness distribution for a uniform disk is simply given by a 2-D top-hat function. 
Thus, the normalized visibilities are given by 

V ( A _ o ^lC^uv) ni v 

* uniform v uv I ^ n i \ ) 

ir6r uv 

where 9 is the angular diameter of the uniform disk brightness distribution, and r uv = 
(u 2 + v 2 ) 1 / 2 is the "uv radius" (Equation 2). The model for the observed squared visibilities 
is obtained by using Equation 6 with V x = Vuniform- 



3.3.3. Accretion Disk Model 



We derive the brightness distribution and predicted visibilities for a geometrically thin 
irradiated accretion disk following the analysis of HSVK and MST. Assuming that the disk 
is heated by stellar radiation and accretion (Lynden-Bell & Pringle 1974), the temperature 
profile (in the regime where R*/R <C 1) is, 



disk 



TlAU (atj) 



-3/4 



(12) 



where T^u is defined as the temperature at 1 AU, given by 



Tiau = 



2-52 xltrM 1^1 T* + 5.27 x 10 10 ( ( , 

RcJ * \M Q J llO-SMoyr- 1 



1/4 



(13) 



We assume that the disk is truncated at an inner radius R in , and an outer radius, i? ou t- 
Guided by Figure 14 of HSVK, we choose i? m to be the radius where the temperature, Tj n , 
is 2000 K. Thus, 

(R \ 3/4 

T 1AU = 2000 . (14) 

2000 K is a likely (upper limit) sublimation temperature for the dust grains that make up 
circumstellar disks, and thus it is reasonable that there would be little or no dust emission 
interior to R m (although the model does not exclude the possibility of optically thin gas 
interior to R in ). We choose R out to be the lesser of 1000 AU or the radius at which T = 3 
K (.Rout is not crucial in this analysis, since most of the near-IR flux comes from the hotter 
inner regions of the disk). 
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The brightness distribution and visibilities for this disk are calculated by determining 
the contributions from a series of annuli from R in to R Q ut- The flux in an annulus specified 
by inner boundary i?j and outer boundary i? f is given by 

Annulus = ^l B M) + B «{Tt)\{R 2 ! ~ (^) 

and the visibilities for this annulus are (following MST), 



Knnulus = -^p[Bv(Ti) + B v (T t )\ 



2 Ji(7r6» f r uv ) 2 Ji(7r6>ir t 



itf It: 



n9 { r uv 1 n9ir v 



(16) 



Here, d is the distance to the source, v is the observed frequency, B u is the Planck function, 
T is the temperature, R is the physical radius, 9 is the angular size, r uv is the "uv radius" 
(Equation 2), and i, f indicate the inner and outer boundaries of the annulus. To obtain the 
visibilities for the entire disk, we sum the visibilities for each annulus, and normalize by the 
total flux: 

ERout t r 
R m V annulus 

Vdisk = _p , ~ • I 1 'J 

^ Rin ^annuius 

The resultant model visibilities are obtained by plugging this expression into Equation 6. 
We note that although we do not use the observed excess K-band flux to constrain the disk 
model, we do verify that the total flux in the model is consistent (to within a factor of 2) 
with the observations. 



3.3.4. Ring Model 
The brightness distribution for a uniform ring model is given by 

I rins (x, y) = { constant if IT < V* 2 + a 2 < ¥ . (18) 
1 otherwise 

Here, (x, y) are the angular offsets from the central star. We define the width of the ring 
via the relation / = W j R, where R is the radius of the inner edge of the ring, and W is the 
width of the ring. Using this relation, we write the inner and outer angular radii of the ring 
as #i n and 6* out = (1 + f)0 in . The normalized visibility of the ring is given by 

King = rre.W + P) [(1 + /} Jl([1 + f ^ 6 ^ ~ J i(^uv)] , (19) 

where r uv = (u 2 + v 2 ) 1 ^ 2 is the "uv radius" (Equation 2). The model for the observed 
visibilities is obtained by using Equation 6 with V x = King- 
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In order to facilitate comparison of our data to puffed up inner disk models from the 
literature, we will use ring widths derived from radiative transfer modeling by DDN. Specif- 
ically, for stars earlier than spectral type B6, we assume / = 0.27, and for stars later than 
B6, we assume / = 0.18 (Table 1 from DDN). 



3.3.5. Two- Component Model 

This model simulates a wide binary, where visibilities are effectively due to two station- 
ary point sources, with some flux ratio and angular separation vector. We explore flux ratios 
from 0.2 to 1, and angular separations from 1 to 100 mas. For flux ratios < 0.2, or angular 
separations < 1 mas, the effects of the companions on the visibilities will be negligible, and 
we can rule out angular separations > 100 mas from adaptive optics imaging (§3.1). The 
squared visibility for the binary model is, 

1 + R 2 + 2.RCOS (^B-s) 
V 2 = V A / (on) 

^binary (1 + i?) 2 ' 

where (B • s)/X — 9[usm(tp) + vcos(ip)], 9 is the angular separation of the binary, ip is the 
position angle, R is the ratio of the fluxes of the two components, and A is the observed 
wavelength. 



3.4. Modeling of Individual Sources 

For each source, we fit the PTI and IOTA visibility data with the models described 
in §3.3 using grids of parameter values. The grid for face-on disk models was generated 
by varying 9 from 0.01 to 10 mas in increments of 0.01 mas. For inclined disk models, in 
addition to varying 9, we varied from 0° to 90° and ip from 0° to 180°, both in increments 
of 1°. As mentioned above, = corresponds to face-on, and ip is measured east of north. 
Since inclined disk models are symmetric under reflections through the origin, we do not 
explore position angles between 180° and 360°. For the binary model, we varied 9 from 1 to 
100 mas in increments of 0.01 mas, ip from 0° to 180° in increments of 1°, and R = F 2 /Fi 
from 0.2 to 1 in increments of 0.001. 

For each point in the parameter grid, we generated a model for the observed u — v 
coverage, and calculated the reduced chi squared (xl) to determine the "best-fit" model. 1- 
<7 confidence limits were determined by finding the grid points where Xr equals the minimum 
value plus one. For inclined disk or binary models, the confidence limits on each parameter 
were determined by projecting the 3-D Xr — min+1 surface onto the 1-D parameter spaces. 
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Tables 4-8 list the best-fit angular size scales (9) for face-on models, the sizes (9), position 
angles (ip), and inclinations (0) for inclined disk models, and the angular separations (9), 
position angles (■0), and brightness ratios (R) for binary models. Values of Xr are a l so 
included in the Tables. Figures 1, 3, 5, 7, and 9 show plots of observed V 2 versus r uv for 
each source along with the curves predicted by various face-on models. Inclined models are 
not circularly symmetric, and the visibilities are a function of the observed position angle 
in addition to the projected baseline (Figure 11). We plot the observed and modeled V 2 for 
inclined models as a function of hour angle in Figures 2, 4, 6, 8, and 10. 

The best-fit binary separations for all sources in our sample are > 2.5 mas. For the 
distances and approximate masses of the sources in our sample, these separations correspond 
to orbital periods of many years. Thus, our assumption that the two point sources in the 
binary model are stationary is justified. 

3.4.1. AB Aur 

The visibilities for AB Aur are consistent with a disk-like circumstellar distribution that 
is inclined by < 30° (Figures 1-2). From Table 4, the best-fit models indicate size scales 4 
between 2.2 and 5.8 mas, and an inclination angle between 27° and 35°. The values of Xr are 
significantly lower for inclined models than for face-on models (x^ ~ 1 an d 2, respectively; 
Table 4), and the data cannot be fit well by a binary model (x 2 r ~ 8). 

3.4.2. VV Ser 

The angular size scales for best-fit disk models range from 1.5 to 3.9, and the disk 
inclinations are between 80° and 90° (Table 5). An inclined disk model clearly fits the VV 
Ser data better than a face-on model (Figures 3 and 4). Inclined model fits give Xr < 1> 
while face-on model fits have xl > 5 (Table 5). However, as indicated in Figure 11, the u — v 
coverage for this object is rather sparse, and precludes placing stringent constraints on the 
value of 0. Moreover, with such sparse u — v coverage a binary model cannot be ruled out 
(Figure 4). 



4 As outlined in §3.3.1-3.3.4, characteristic size scales for different models measure different parts of the 
brightness distributions: Gaussian models measure full widths at half maxima, uniform disk models measure 
outer diameters, accretion disk models measure inner disk diameters, and ring models measure inner ring 
diameters. The spread in quoted angular sizes for a source is mainly due to these differences. 
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3.4-3. VI 68 5 Cyg 

The size scales for V1685 Cyg under the assumptions of various disk models range from 
1.3 to 3.9 mas, and the inclinations are between 49° and 51° (Table 6). The visibility data are 
not fit very well by any model, although of those considered, inclined disks fit best (Figures 
5-6). While we cannot rule out a binary model, we note that xl ~ 3 for the binary model, 
compared to xl ~ 2 for inclined disk models (Table 6). Better coverage of the u — v plane 
should help to improve our understanding of this source (Figure 11). 

3.44. AS 442 

The PTI data for AS 442 generally have low signal-to-noise, and it is difficult to distin- 
guish between different models. Nevertheless, we can make an approximate determination of 
the size scale, although we cannot distinguish between inclined disk, face-on disk, or binary 
models (Figures 7 and 8). The size scales for various disk models range from 0.9 to 2.7 mas 
(Table 7). 

3.4.5. MWC 1080 

The PTI and IOTA observations for MWC 1080 are completely incompatible with face- 
on models (Xr > 40), and significantly non-zero inclinations are required to fit the data well 
(Figures 9 and 10). The best-fit inclination angles for various disk models range from 51° 
to 56°, and the angular size scales are between 1.5 and 4.1 mas (Table 8). For this source, 
we can rule out a binary model with a relatively high degree of confidence: Xr ~ 10 f° r the 
binary model, compared to Xr ~ 2 for inclined disk models. 

4. Discussion 

As discussed in §1, there is currently a wide variety of evidence that supports the 
existence of circumstellar disks around many HAEBE stars. Our new PTI results strengthen 
this contention. Resolved, small-scale (~ 1 AU) distributions of dust are found in all observed 
sources, and the non-symmetric intensity distributions of best-fit models for most objects 
provide support for inclined disk geometries. 

We suggest that the material around VV Ser, V1685 Cyg, and MWC 1080 is significantly 
inclined, and we cannot rule out a high inclination angle for AS 442. This hypothesis is com- 
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patible with observed optical variability in VV Ser and AS 442 (AVyvser ~ 2, AVas442 ~ 1; 
Herbst & Shevchenko 1999), which has been attributed to variable obscuration from clumps 
of dust orbiting in inclined circumstellar disks. The AB Aur data, in contrast, are consistent 
with a circumstellar distribution that is within 35° of face-on. This agrees well with MST 
and is compatible with modeling of scattered light observed with the Hubble Space Tele- 
scope, which suggests an inclination angle < 45° (Grady et al. 1999). The small amplitude 
of variability in AB Aur (Ay ~ 0.25; Herbst & Shevchenko 1999) is also consistent with this 
low inclination angle (under the assumption that variability is caused by time-dependent cir- 
cumstellar obscuration). The low inclination angle does not, however, agree with mm- wave 
imaging in the 13 CO (1-0) line, which yields an estimated inclination of 76° for the AB Aur 
disk (Mannings & Sargent 1997). 

The angular sizes determined from our observations are generally in good agreement 
with the non-inclined (0 = 0) flat accretion disk models of HSVK for early-type Herbig Be 
stars, V1685 Cyg and MWC 1080, but not for the later-type stars, AB Aur, VV Ser, and 
AS 442 (the spectral type for VV Ser, AO, is uncertain by ±5 spectral subclasses; Mora et 
al. 2001). Angular sizes derived from the earlier IOTA observations (MST) were often an 
order of magnitude larger than those predicted by the HSVK models, and on this basis MST 
ruled these models out. 

HSVK determined the best-fit models for the SEDs of HAEBE sources by assuming 
a face-on disk geometry, adjusting the accretion rate to match the mid-IR flux, and then 
adjusting the size of the inner hole to match the near-IR flux. We compare our results 
with theirs in a qualitative way by plotting the visibilities predicted by the HSVK models 
along with the observed PTI and IOTA visibilities in Figures 1, 3, 5, 7, and 9. For a 
more quantitative comparison, we use published luminosities, effective temperatures, and 
accretion rates (HSVK) to calculate the inner radii predicted by flat accretion disk models 
with T in = 2000 K (Equations 13 and 14), and compare these estimates to our interferometric 
results (which were also derived assuming T in = 2000 K; §3.3.3). In Table 9, -Rf aC c-on and 
-Rinciined are the inner radii determined by fitting the interferometric data to face-on and 
inclined accretion disk models, respectively (§3.3.3), and Rm =0 , -Rm^o are ^ ne ra dii calculated 
using the HSVK flat disk models without accretion, and with accretion effects included, 
respectively. No estimate of M is available for AS 442. 

Our data for the later-type stars AB Aur, VV Ser, and AS 442 are fairly consistent with 
the puffed up inner disk models of DDN, assuming inner disk temperatures > 2000 K. In 
contrast, puffed up inner disk models are completely incompatible with the PTI results for 
the very early-type stars in our sample, V1685 Cyg and MWC 1080. The radius of the inner 
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wall, Ri n , predicted by DDN is, 

R * = {Sp+fi' < 21) 

where, L* is the (published) stellar luminosity, T in is the temperature of the inner wall, and / 
is the ratio of the width of the inner wall to its radius. Based on DDN, we assume / = 0.27 
for stars earlier than spectral type B6, and / = 0.18 for later-type stars. We calculate 
i?i n for Ti n = 1500, 2000 K (likely sublimation temperatures for silicate and graphite dust 
grains, respectively), and compare these to the ring diameters derived from fitting to near- 
IR interferometric visibilities. In Table 10, -Rf aC c-on and inclined represent the inner radii 
determined for face-on and inclined ring models, respectively (§3.3.4), and -fi^ootb -R1500 are 
the inner radii predicted by the DDN puffed-up inner disk models, assuming sublimation 
temperatures of 2000 and 1500 K, respectively. 

We note that the comparison of our interferometric results to physical models should be 
independent of the assumed distance (see Appendix A). The inner radius is oc L 1 ^ 2 oc d in 
both the DDN and HSVK models, and the linear sizes determined from our interferometric 
results (converted from modeled angular sizes) are also oc d, and thus, the comparison is 
independent of d. 

Flat accretion disk models (HSVK) are generally in good agreement with the observed 
visibility data for early-type B-stars, while puffed up inner disk models (DDN) seem more 
consistent for later-type stars. We speculate that this could be due to different accretion 
mechanisms in earlier and later-type stars. A similar idea has been put forward based on the 
results of Ha spectropolarimetry, where differences in the observations for early-type HBc 
stars and later-type HAe stars have been attributed to a transition from disk accretion in 
higher- mass stars to magnetic accretion in lower-mass stars (Vink et al. 2003). 

There is always the possibility that the visibilities for some of the observed HAEBE 
sources may be (partially) due to close companions. For AB Aur and MWC 1080, we can 
rule out binary models (with separations > 1 mas) with a high degree of confidence. However, 
MWC 1080 is an eclipsing binary with a period of P tv 2.9 days (Shevchencko et al. 1994; 
Corporon & Lagrange 1999). The separation is much too small to be detected by PTI, and 
the observed visibilities for this source are thus probably due to an inclined circum- binary 
disk. Observations over a time-span of ~ 100 days (Table 2), with visibilities that are fairly 
constant in time (Figures 5 and 6) provide some evidence against V1685 Cyg being a binary. 
As yet, the binarity status of AS 442 and VV Ser remain uncertain based on our visibility 
data, although radial velocity variations of spectral lines in AS 442 have been attributed to 
a binary with P 64 days and e 0.2 (Corporon & Lagrange 1999). 
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5. Summary 

We observed the HAEBE sources AB Aur, VV Ser, V1685 Cyg (BD+40°4124), AS 442, 
and MWC 1080 at 2.2 /zm with the Palomar Testbed Interferometer. These are only the 
second published near-IR interferometric observations of HAEBE stars. From these high 
angular resolution data, we determined the angular size scales and orientations predicted 
by uniform disk, Gaussian, ring, and accretion disk models, and we examined whether the 
data were consistent with binary models. AB Aur appears to be surrounded by a disk that 
is inclined by < 30°, while VV Ser, V1685 Cyg, and MWC 1080 are associated with more 
highly inclined circumstellar disks. With the available data, we cannot distinguish between 
different radial distributions, such as Gaussians, uniform disks, rings, or accretion disks. 

While the angular size scales determined in this work are generally consistent with the 
only other near-IR interferometric measurements of HAEBE stars by MST, our measure- 
ments are the first that show evidence for significantly inclined morphologies. Moreover, the 
derived angular sizes for early type Herbig Be stars in our sample, V1685 Cyg and MWC 
1080, agree fairly well with those predicted by face-on accretion disk models used by HSVK 
to explain observed spectral energy distributions. The observations of AB Aur, VV Ser, 
and AS 442 are, however, not entirely compatible with these models, and may be better 
explained through the puffed-up inner disk models of DDN. 
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A. Distance Estimates 



AB Aur is associated with the Taurus-Auriga molecular cloud, and thus the estimated 
distance to this source (d = 140 pc) is accurate to ~ 10%. Photometric studies of VV 
Ser and other stars in Serpens estimated distances of d ~ 250 pc (Chavarria-K. et al. 1988) 
and d « 310 pc (de Lara et al. 1991), while an earlier study based on photometry of a 
single source estimated d = 440 pc (Strom et al. 1974). Based on these estimates, we adopt 
a distance of 310 pc. Distance estimates to V1685 Cyg range from 980 pc (based on an 
extinction-distance diagram for 132 stars within 3.5°; Shevchencko et al. 1992) to 1000 pc 
(based on locating V1685 Cyg on the main sequence; Strom et al. 1972), to 2200 pc (based 
on photometry of stars in a large scale region around V1685 Cyg; Hiltner & Johnson 1956). 
We adopt a distance of d = 1000 pc to V1685 Cyg, since the 2200 pc estimate would imply a 
luminosity higher than expected for the published spectral type. AS 442 is associated with 
the North American Nebula, and thus the adopted distance of 600 pc is probably accurate to 
~ 10%. The distance to MWC 1080 has been determined by fitting photometric observations 
to the main sequence (d = 1000 pc; HSVK), and using the Galactic rotation curve (d = 2500 
pc; Canto et al. 1984). We adopt a distance of 1000 pc to MWC 1080, since the 2500 pc 
estimate based on the Galactic rotation curve would imply a luminosity much higher than 
expected for the published spectral type. Moreover, the 2500 pc estimate is uncertain by 
~ 50%, while the 1000 pc estimate is accurate to ~ 20%. 
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Fig. I.— V 2 data from PTI (symbols) and IOTA (filled dots; MST) for AB Aur, as a 
function of r uv = [u 2 + v 2 ) 1 ^ 2 . PTI data for individual nights are represented by different 
symbols. Face-on uniform disk (solid line), ring (dotted line), and accretion disk (dashed 
line) models are over-plotted. We also plot the visibilities calculated for an accretion disk 
model with R m = 0.09 AU and T in = 2360 K (HSVK; thick solid line). 
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Fig. 2. — PTI and IOTA V 2 data for AB Aur (represented as in Figure 1), as a function of 
hour angle. Over-plotted are face-on and inclined accretion disk models (solid and dotted 
lines, respectively), as well as the best-fit binary model (dashed line). 
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Fig. 3. — PTI V 2 data for VV Ser, as a function of r uv = (-u 2 + t> 2 ) 1//2 . PTI data for individual 
nights are represented by different symbols. Face-on uniform disk (solid line), ring (dotted 
line), and accretion disk (dashed line) models are over-plotted. We also plot the visibilities 
calculated for an accretion disk model with R in = 0.08 AU and T in = 2710 K (HSVK; thick 
solid line). 
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Fig. 4. — PTI V 2 data for VV Ser (represented as in Figure 3), as a function of hour angle. 
For clarity, we have plotted V 2 + 1 for the data taken with the NW baseline. Over-plotted 
are face-on and inclined accretion disk models (solid and dotted lines, respectively), as well 
as the best-fit binary model (dashed line). 
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Fig. 5.— PTI V 2 data for V1685 Cyg, as a function of r uv = {u 2 + v 2 ) 1 ' 2 . PTI data for 
individual nights are represented by different symbols. Face-on uniform disk (solid line), 
ring (dotted line), and accretion disk (dashed line) models are over-plotted. We also plot 
the visibilities calculated for an accretion disk model with R in = 0.63 AU and T in = 2060 K 
(HSVK; thick solid line). 
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Fig. 6. — PTI V 2 data for V1685 Cyg (represented as in Figure 5), as a function of hour 
angle. For clarity, we have plotted V 2 + 1 for the data taken with the NW baseline. Over- 
plotted are face-on and inclined accretion disk models (solid and dotted lines, respectively), 
as well as the best-fit binary model (dashed line). 
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Fig. 7. — PTI V 2 data for AS 442, as a function of r uv 
nights are represented by different symbols. Face-on uniform disk (solid line), ring (dotted 
line), and accretion disk (dashed line) models are over-plotted. We also plot the visibilities 
calculated for an accretion disk model with R 1D = 0.10 AU and T in = 2000 K (HSVK; thick 
solid line). 
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Fig. 8. — PTI V 2 data for AS 442 (represented as in Figure 7), as a function of hour angle. 
For clarity, we have plotted V 2 + 1 for the data taken with the NW baseline. Over-plotted 
are face-on and inclined accretion disk models (solid and dotted lines, respectively), as well 
as the best-fit binary model (dashed line). 
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Fig. 9.— V 2 data from PTI (symbols) and IOTA (filled dots; MST) for MWC 1080, as a 
function of r uv = (u 2 + v 2 ) 1 ^ 2 . PTI data for individual nights are represented by different 
symbols. Face-on uniform disk (solid line), ring (dotted line), and accretion disk (dashed 
line) models are over-plotted. We also plot the visibilities calculated for an accretion disk 
model with R in = 0.59 AU and T in = 2490 K (HSVK; thick solid line). While none of these 
face-on models fit the data well, good fits are obtained with inclined models (Figure 10). 
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Fig. 10.— PTI and IOTA V 2 data for MWC 1080 (represented as in Figure 9), as a function 
of hour angle. Over-plotted are face-on and inclined accretion disk models (solid and dotted 
lines, respectively), as well as the best-fit binary model (dashed line). Note the significant 
improvement in the fit when inclination effects are included in the model. 
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Fig. 11. — Contour plots of best-fit inclined uniform disk models for AB Aur, VV Ser, V1685 
Cyg, AS 442, and MWC 1080, whose parameters are listed in Tables 4-8. The contour 
increment is 10% in V 2 . We also plot the best-fit inclined disk model for MWC 1080 on 
the sky (bottom right panel, greyscale). We over-plot the uv points sampled for each source 
by the PTI NW baseline (open triangles), the PTI NS baseline (open diamonds), and by 
IOTA (filled dots). Since we know that the brightness distributions of the sources are real, 
the visibilities must be reflection symmetric (through the origin), and so we also plot these 
reflections of the sampled uv points. 



Table 1. Observed Sources 



Source 


Alt. Name 


a (J2000) 


5 (J2000) 


d (pc) 


Sp.Ty. 


V 


K 


Ft (Jy) 


Fl (Jy) 


AB Aur 


HD 31293 


04 h 55 m 45.84 s 


+30°33'04"3 


140 


AOpe 


7.07 


4.27 


1.92 


10.59 


VV Ser 


HBC 282 


18 h 28 m 49.00 s 


+00°08'39"0 


310 


AOVevp 


11.90 


6.44 


0.20 


1.85 


V1685 Cyg 


BD+40°4124 


20 h 20 m 28.25 s 


+41°21'51"6 


1000 


B2Ve 


10.71 


5.70 


0.42 


3.64 


AS 442 


V1977 Cyg 


20 h 47 m 37.47 s 


+43°47'24"9 


600 


B8Ve 


10.89 


6.75 


0.20 


1.21 


MWC 1080 


V628 Cas 


23 h 17 m 26.10 s 


+60°50'43"0 


1000 


BOeq 


11.68 


4.83 


0.87 


9.85 



References. — Distances, spectral types and V magnitudes from Hillenbrand et al. (1992), Mora et al. (2001), 
Strom et al. (1972), de Lara et al. (1991), and Bigay & Garner (1970). For discussion of the adopted distances, 
see Appendix A. K magnitudes from the present work, f : De-reddened fluxes. 



- 32 - 





Table 2. 


Summary of Observations 




OU 111 L^t; 


YJ d b L, I 1V1 YJ I 


Baseline 


ha cov." 1 " 


vjdll UI dlUl o 1 YYYJ J 


AB Aur 


52575 


NW 


[1.21,1.85] 


29645, 32301 




52602 


NW 


[-1.95,1.51] 


29645, 32301 


VV Ser 


52490 


NW 


[-0.72,0.54] 


171834 




52491 


NS 


[-1.55,-0.74] 


171834 




52493 


NW 


[-1.31,0.20] 


171834 




52499 


NW 


[-0.96,0.84] 


164259,171834 


V1685 Cyg 


52418 


NW 


[-1.10,-1.00] 


192640,192985 




52475 


NW 


[-1.69,1.44] 


192640,192985 




52476 


NW 


[-1.80,-0.48] 


192640,192985 




52490 


NW 


[-0.97,1.70] 


192640,192985 




52491 


NS 


[-1.27,2.38] 


192640 




52492 


NW 


[-0.90,-0.90] 


192640 




52545 


NS 


[-1.12,2.48] 


192640,192985 


AS 442 


52475 


NW 


[0.21,1.33] 


192640,192985 




52476 


NW 


[-0.21,-0.21] 


192640,192985 




52490 


NW 


[-1.11,0.38] 


192640 




52491 


NS 


[-0.69,2.54] 


192640 




52492 


NW 


[-1.05,0.00] 


192640 




52545 


NS 


[-0.87,1.58] 


192640,192985 


MWC 1080 


52475 


NW 


[0.17,0.17] 


219623 




52476 


NW 


[-1.99,0.52] 


219623 




52490 


NW 


[-0.14,1.39] 


219623 



References. — f: Hour angle coverage of the observations. 
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Table 3. Properties of Calibrator Sources 



Name 


a (J2000) 


5 (J2000) 


Sp.Ty. 


V 


K 


Cal. Size (mas) 


Aa n 


HD 29645 


04 h 41 m 50.26 s 


+38°16'48"7 


GOV 


6.0 


4.6 


0.56 ±0.09 


8.2 


HD 32301 


05 h 03 m 05.75 s 


±21°35'23"9 


A7V 


4.6 


4.1 


0.47 ±0.10 


9.1 


HD 164259 


18 h 00 m 29.01 s 


-03°41'25"0 


F2IV 


4.6 


3.7 


0.77 ±0.08 


7.5 


HD 171834 


18 h 36 m 39.08 s 


+06°40'18"5 


F3V 


5.4 


4.5 


0.54 ±0.07 


6.8 


HD 192640 


20 h 14 m 32.03 s 


+36°48'22"7 


A2V 


4.9 


4.9 


0.46 ±0.02 


4.7\9.4 2 


HD 192985 


20 h 16 m 00.62 s 


+45°34'46"3 


F5V 


5.9 


4.8 


0.44 ±0.04 


4.3\5.9 2 


HD 219623 


23 h 16 m 42.30 s 


+53°12'48"5 


F7V 


5.6 


4.3 


0.54 ±0.03 


9.5 



References. - 



1,2: Offsets from V1685 Cyg, AS 442, respectively. 
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Table 4. Results of Modeling for AB Aur 



Model Xr (mas) (°) (°)t 



Face-On Gaussian 


2.42 


O.O»_ Qg 






Face-On Uniform 


1.78 


r o 4 +0.09 






Face-On Accretion 


1.93 


o I o+0.03 
z - io -0.02 






Face- On Ring 


1.92 


o 9fi+ - 02 

°- z,{J -0.03 






Inclined Gaussian 


0.96 


q OO+0.38 


103 + | 


35111 


Inclined Uniform 


0.89 


5.801851 


128+3° 


2611° 


Inclined Accretion 


0.88 


2-301^ 


10512$ 


07+13 
z '-17 


Inclined Ring 


0.88 


°- uu -0.38 


1441^ 


2811° 


Binary Model 


8.96 


o 41 +0.13 
°-^ ± -0.28 


38+J 


0-581°;^ 



References. - f: For the binary model, represents 
the brightness ratio, R = F 2 / ' F\. 
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Table 5. Results of Modeling for VV Ser 



Model Xr (mas) (°) (°)t 



Face-On Gaussian 


9.13 


o qq+0.09 
z -°°-0.09 






Face-On Uniform 


6.91 


3.6818:S 






Face-On Accretion 


8.33 


1-49183S 






Face- On Ring 


5.86 


9 qn+0.07 
z - ou -0.07 






Inclined Gaussian 


0.85 


9 ^f, +lm 
z - du -0.13 


37lJ 5 


891^ 


Inclined Uniform 


0.85 


3.94+2-35 


41 ;f 


82ll 3 


Inclined Accretion 


0.85 


1 62 +L58 

-■-• uz -0.98 


38_ 70 


83-45 


Inclined Ring 


0.85 


2.441J;?? 






Binary Model 


0.85 


8.8oiS:§§ 


17611 


0.4518:82 



References. - f : For the binary model, represents 
the brightness ratio, R = F 2 /Fi. 
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Table 6. Results of Modeling for V1685 Cyg 



Model Xr (mas) (°) (°)t 



Face-On Gaussian 


6.51 


1 Q6+ 011 
i - yo -o.n 






Face-On Uniform 


7.52 


q i 7+O.I6 
°- i ' -0.15 






Face-On Accretion 


6.75 


i 9 7+0.06 
L - z ' -0.06 






Face- On Ring 


8.10 


i Q2+ - 09 






Inclined Gaussian 


2.32 


9 a O+0.44 
z -+°-0.37 


1257^28 


5H12 


Inclined Uniform 


2.36 


o qi+o.60 

°- y± -0.55 


124^ 4 


50l}i 


Inclined Accretion 


2.33 


i cr 7 +0.27 
l .O 1 _ 22 






Inclined Ring 


2.38 


9 qq+0.37 
Z.OO_ 29 


122152 


491S 


Binary Model 


3.33 


o 4 i+0.37 
°-^ i -0.56 


62^ 


0.2418:82 



References. - f: For the binary model, represents 
the brightness ratio, R = F 2 / F\. 
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Table 7. Results of Modeling for AS 442 



Model Xr (mas) ip (°) (°)t 



Face-On Gaussian 


0.99 


1 4O+0-19 
i -^ y -0.19 






Face-On Uniform 


1.04 


2 44+ - 29 

z -^-0.28 






Face-On Accretion 


0.99 


0.95183 






Face-On Ring 


1.07 








Inclined Gaussian 


0.94 


l.UO_ 29 


60 +12 ° 

uu -60 


4i -41 


Inclined Uniform 


0.94 


2 R7+ 1 - 29 
z -° ' -0.34 


°°-63 


oy -39 


Inclined Accretion 


0.94 


i - uo -0.18 


°°-63 


°°-36 


Inclined Ring 


0.95 


1.70lg;g 


fi c+115 


qo+52 


Binary Model 


0.95 


2.691?;" 


30lg 


21+ 79 
u - z± -o.oi 



References. - f: For the binary model, represents 
the brightness ratio, R — F 2 / F±. 
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Table 8. Results of Modeling for MWC 1080 



Model Xr (mas) ip (°) (°)t 



Face-On Gaussian 


56.33 


9 O/1+0.05 
z -°^-0.05 






Face- On Uniform 


42.04 


3.8418$ 






Face-On Accretion 


54.24 


i-54l8:S 






Face-On Ring 


36.00 


9 qo+0.04 
z -°°-0.05 






Inclined Gaussian 


3.21 


2 61+ 011 

^• ul -0.08 


n +u 

' L -9 


56+6 


Inclined Uniform 


2.54 


4 1 s+ - 12 

4 - i<3 -0.10 






Inclined Accretion 


3.07 


1 6Q+ - 07 

i - oy -0.05 


71 + 10 


55+1 


Inclined Ring 


2.28 


9 /| 7+O.O6 
-0.06 




51±S 


Binary Model 


9.32 


9 p-7+0.22 
z - d ' -0.18 


56* 


0.36±g;g 



References. — f: For the binary model, represents the 
brightness ratio, R — F 2 /Fi. 
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Table 9. Comparison with Hillenbrand et al. (2001) Models 



Source 


p* 

■^f ace— on 


Ft* 

-^inclined 








(AU) 


(AU) 


(AU) 


(AU) 


AB Aur 


0.15 ±0.01 


0.16 ±0.01 


0.07 


0.12 


VV Ser 


0.23 ±0.01 


0.25 ±0.19 


0.03 


0.13 


V1685 Cyg 


0.64 ±0.03 


0.79 ±0.13 


0.44 


0.71 


AS 442 


0.29 ±0.04 


0.31 ±0.12 


0.10 




MWC 1080 


0.77 ±0.02 


0.85 ± 0.03 


0.79 


0.79 



References. — *: Error bars based on 1-a uncertainties 
of best-fit face-on and inclined accretion disk models. 
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Table 10. Comparison with Dullemond et al. (2001) Models 



Source 


face— on 


p* 

inclined 


-R2000 


-R1500 




(AU) 


(AU) 


(AU) 


(AU) 


AB Aur 


0.23 ±0.01 


0.26 ±0.03 


0.32 


0.57 


VV Ser 


0.36 ±0.01 


0.38 ±0.14 


0.24 


0.42 


V1685 Cyg 


0.96 ±0.05 


1.17 ± 0.16 


3.71 


6.59 


AS 442 


0.47 ±0.05 


0.56 ±0.13 


0.52 


0.93 


MWC 1080 


1.17±0.02 


1.24 ±0.03 


8.69 


15.45 



References. — *: Error bars based on 1-a uncertainties 
of best-fit face-on and inclined ring models. 



